Double Beta Decay of 100 Mo to Excited Final States 



M. J. HornislE and L. De BraeckeleeiQ 

The Department of Physics, 
Duke University and Triangle Universities Nuclear Laboratory, 
Durham, NC 27708-0308 



in 
o 
o 

(N 
O 

Q 

o 

(N 



A. S. Barabash and V. I. Umatov 
Institute for Theoretical and Experimental Physics, 
B.Cheremushkinskaya 25, 117259, Moscow, Russia 
(Dated: February 8, 2008) 

A systematic study of the inclusive (Ou + 1v) double beta (/3/3) decay of 100 Mo to various ex- 
cited final states of 100 Ru was performed. Utilizing two large HPGe detectors operated in coin- 
cidence, a search for the subsequent deexcitation 77 cascades was conducted. A 1.05-kg sample 
of isotopically enriched (98.4%) 100 Mo was investigated for 455 days, yielding an unambiguous 
observation of the /3/3 decay of 100 Mo to the Of state (1130.3 keV) of 100 Ru. This excited final 
state decays via the 0^ — > 2^ — > 0^ s sequence (with E^x = 590.8 keV and E l2 = 539.5 keV), 
and 22 such coincidence events were detected, with a continuous background estimated to be 2.5 
events. This counting rate corresponds to a decay half-life for the /3(3(0 + — > 0^~) transition of 
T^2 + = [6.0l^'j(stat) ± 0.6(syst)] x 10 20 years. Lower limits on decay half-lives were achieved 
for higher excited final states. 

PACS numbers: 21.10.Tg, 23.40.-s, 27.60. +j 
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I. INTRODUCTION 

A considerable escalation in the attention paid to the 
study of double beta (/?/?) decay is owed, in no small 
part, to the tremendous progress made by recent neu- 
trino oscillation experiments. The remarkable success of 
solar and atmospheric neutrino experiments (see reviews 
represents a profound achievement and a verita- 
ble breakthrough in our understanding of the neutrino. 
More recently, the KamLAND reactor antineutrino result 
and the latest SNO result [j| further reinforce the ev- 
idence for the oscillation of neutrinos. Indeed, there can 
no longer be any doubt that the neutrino is, in fact, a 
massive particle. 

At the same time, however, for all of the progress that 
has been accomplished in measuring the mass differences 
(Am 2 ) for neutrinos of unlike flavors, such oscillation ex- 
periments are inevitably and inherently insensitive to two 
additional important properties: (1) the neutrino's abso- 
lute mass scale and (2) the Majorana-Dirac nature of the 
neutrino. The only present experimental method capa- 
ble of attempting to simultaneously answer both of these 
questions is the search for the neutrinoless mode (Ou/3/3) 
of double beta decay. Although this transition has not yet 
been observed beyond a reasonable doubt, the successful 
detection of this decay mode may also provide useful in- 
formation about the neutrino mass hierarchy (normal, 
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inverted, or quasi-degenerate) and lepton-sector CP vio- 
lation via a measurement of the Majorana CP-violating 
phases (see discussions 0,0,13). 

The recent surge in interest in j3[3 decay has ini- 
tiated the proposal of several large-mass experiments 
that will search for the neutrinoless mode (see review 
[HI EH H3 ) ^he interpretation of any future results 
from these experiments depends on the reliable estima- 
tion of nuclear matrix elements (M 0v ) for the OisfiP tran- 
sition. While the phase-space factor for this decay can be 
accurately calculated, any uncertainty in the correspond- 
ing nuclear matrix elements will adversely affect the abil- 
ity to correctly extract Majorana neutrino mass informa- 
tion. Several theoretical models have been devised, in- 
cluding the well-known quasiparticle random phase ap- 
proximation (QRPA), but a thorough analysis of the ac- 
curacy of the various models to properly calculate the 
nuclear matrix elements is difficult at present, although 
efforts are being made (see, e.g., fl3lll4) '). 

One testing ground for our understanding of the /3/3- 
decay nuclear matrix elements is through the study of 
the ordinary allowed second-order weak decay, namely 
the two-neutrino double beta (2v(3(3) decay. In general, 
expanding our experimental knowledge of the 2v(3j3 pro- 
cess for a variety of different nuclei will improve the over- 
all comprehension of the nuclear part of double beta de- 
cay. Although the nuclear matrix elements for this mode 
{M 2v ) are not identical to those of the neutrinoless mode, 
experimental measurements of M 2u can be compared to 
theoretical model calculations as a check of the validity 
of the various theoretical schemes upon which the models 
rely. 

Over the past decade and a half, significant progress 
has been made in successfully detecting the 2vf}j3 de- 
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FIG. 1: Simplified /3/3-decay scheme of 100 Mo to the 0+ ex- 
cited state, whose deexcitation proceeds 100% of the time 
through the intermediate 2~f state. 



cay to ground states for several nuclei, including 8 Ca, 
76 Ge, 100 Mo and 150 Nd (see summary [15J). However, 
for certain nuclei, there are low-lying excited levels in the 
daughter nuclide that may also be populated by 0/3 de- 
cay. Although these excited-state transitions suffer from 
a reduced phase space, and thus a longer half-life, it is 
important to realize that such decays present an addi- 
tional quantity to measure, namely the deexcitation 7 
rays. In fact, in the framework of QRPA models, excited- 
state transitions exhibit a much different dependence on 
the particle-particle strength parameter g pp than ground- 
state transitions |16l Il7l Il8| . Hence, different aspects of 
the theoretical techniques used to calculate nuclear ma- 
trix elements may be probed by the study of decays to 
excited states. 

For the case of 100 Mo, the Q values for the /3/3 decay 
to some of the excited states of 100 Ru are large enough 
that the detection of the 2v{3(3 decay may be possible. 
The idea 0] to search for the 7 rays emitted from the 
deexcitation of the excited daughter nucleus prompted 
several measurements, which used enriched molybdenum 
samples, to operate 7-ray detectors in singles mode in 
underground, very-low background settings. By detect- 
ing the subsequent gamma decay of the excited daughter 
nuclide rather than the two beta particles, this type of 
measurement is not sensitive to the particular mode of 
decay. However, the emitted 7 rays from the excited 
nucleus have a fixed energy, so if they can be detected 
with good energy resolution and relatively high efficiency, 
the measurement is rather straightforward and the rele- 
vant experimental background can be substantially sup- 
pressed. 

In particular, the study of 100 Mo focused primarily on 
the 2v(3j3 decay to the 0+ state of 100 Ru (1130.3 keV) 



by searching for the subsequent 0^ — > 2f — > 0+ s decay 
cascade, with E y i = 590.79 keV and E l2 = 539.51 keV 
(see Figure nj. One such experiment was performed by 
measuring enriched 100 Mo samples (~ 1 kg) with a low- 
background HPGe detector. An analysis of the single 
7-ray spectrum positively identified the 2v(3j3(0 + — > Of) 
transition, and a deduced half- life of T\/ 2 = (6.lj^"^) x 
10 20 years was produced [20j. Afterwards, a second ex- 
periment using the same technique yielded a measured 
half-life for this transition of T 1/2 = (9.3±?;f) x 10 20 
yea rs, with a systematic error estimated to be ~ 15% 
|2l| . By summing the single 7-ray spectrum of Ref. [2(j 
to that of Ref. 21], one obtains a combined half- life of 
^1/2 = (7.6 jv? ) x 10 20 years, also with a systematic er- 
ror of 15% |2l|. On a related note, it is interesting to 
point out that a recent experiment using a similar tech- 
nique successfully observed the 2v(3(3 decay of 150 Nd to 
the 0^~ state of 150 Sm, with a half-life estimated to be 
T 1/2 = [l.4±g:|(stat) ± 0.3(syst)] x 10 20 years j2|. 

The preliminary results from the present experiment 
on 100 Mo, with a quoted half-life of this 2i//3/3(0+ -> Of) 
transition of T 1/2 = [5.9±^(stat) ± 0.6(syst)] x 10 20 
years, were previously presented |23j and will be explored 
in greater detail in this Article. 



II. EXPERIMENTAL DETAILS 

By relying on single 7-ray searches, previous mea- 
surements of the 2v(3ft{0 + — > Of ) decay of 100 Mo have 
involved the implementation of low-background detec- 
tion systems, where detectors were made from low- 
radioactivity materials. These experiments were oper- 
ated in underground laboratories, where the natural over- 
head shielding reduces cosmic-ray induced backgrounds. 
These experiments produce single 7-ray spectra exhibit- 
ing an excess of counts in the regions of interest near 
539.5 keV and 590.8 keV, above a rather substantial back- 
ground. In the two cases mentioned in Section [I] the 
correspond ing signal-to-background ratios were approxi- 
mately 1:7 |2£| and 1:4 [2]]]. Extracting the total number 
of counts attributed to the f3(3 decay depends greatly on 
a reliable fit to the background. 

The technique utilized in the present Article involves 
the use of two low-background HPGe detectors operated 
in a coincidence scheme, the first time such a 77 coinci- 
dence technique has been applied to the study of (3f3 de- 
cay. In some cases, as in the present one, the reduced de- 
tection efficiency of the coincidence measurement can be 
more than compensated by a corresponding suppression 
of the associated background. This alternative approach 
to background reduction relies on the simultaneous de- 
tection of the two 7 rays emitted from the deexcitation 
of the daughter nuclide. Owing to the unique 7-ray en- 
ergies, the coincidence background can be suppressed to 
such an extent that the signal-to-background ratio may 
be much higher than in single 7-ray searches, even in 
above-ground experiments. 
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A. TUNL-ITEP Apparatus 

Conducted in the Low Background Counting Facility 
of the Triangle Universities Nuclear Laboratory (TUNL), 
the experimental work is centered around a 77 coinci- 
dence apparatus, the cornerstone of which is the op- 
eration in coincidence of two high-purity germanium 
(HPGe) detectors. These large, custom-made HPGe de- 
tectors (p-type coaxial) were specially designed to have a 
large frontal surface area. The size of each HPGe crystal 
was 8.8 cm in diameter and 5.0 cm in thickness. Each 
detector was fabricated from low-background materials, 
and each crystal was coupled to a very-low background 
cryostat in the J-type configuration, which is used not 
only to isolate the detector crystal from the preamplifier, 
the HV filter and LN2 dewar, but also to allow for its 
easy insertion into an anticoincidence annulus. The two 
HPGe detectors have 1.8-keV FWHM energy resolution 
at 1.33 MeV and 0.8-keV resolution at 0.122 MeV, with 
the relevant peak-to-Compton ratio rated as 80.4. 

To suppress the inherent background contributions 
that result from primary cosmic-ray radiation, secondary 
cosmogenically-induced radiation and primordial radia- 
tion, active and passive shielding techniques were em- 
ployed. The active shielding consists of anticoincidence 
counters, which include two plastic plate scintillators (10- 
cm thickness) and one Nal(Tl) annulus (56-cm length, 
35.6-cm diameter with a 12.5-cm hole along the axis of 
symmetry) into which the source and HPGe detectors are 
inserted. These counters also veto potential backgrounds 
arising from unwanted Compton scattering events within 
the HPGe detectors. The passive shielding includes a 
lead-brick enclosure and the overhead shielding afforded 
by the room ceiling and the floors of the building (~10 
m w.e.). A schematic view of the apparatus can be seen 
in Figure [2 

A reliable energy calibration of the HPGe detectors 
was attained by monitoring known background 7 peaks 
[238.6 keV ( 212 Pb), 511.0 keV (annihilation), 609.3 keV 
( 214 Bi), and 1460.8 keV ( 40 K)]. Various performance 
checks were used to verify the proper operation of the 
detectors and the electronics over the course of the ex- 
periment. For example, daily monitoring of the location 
of the more prominent background 7 peaks allowed for 
any corrections due to gain shifts. Secondly, the sensi- 
tivity of the veto system was readily checked by mon- 
itoring the strength of known single gamma lines (e.g., 
1461-keV 7 rays from 40 K), which in principle are not 
vetoed, against those emitted in coincidence with addi- 
tional gamma quanta (e.g., 511-keV annihilation 7 rays), 
which will normally be vetoed. 

For the study of the 2vf3j3 decay of 100 Mo to excited 
states of 100 Ru, an enriched metallic disk was fabri- 
cated and loaned to TUNL by Russian collaborators from 
the Institute of Theoretical and Experimental Physics 
(ITEP). It was enriched in 100 Mo to 98.4% and had a 
mass of 1.05 kg. The disk was approximately 1.1 cen- 
timeters thick and 10.6 centimeters in diameter, which 




FIG. 2: (Color online) Schematic view (not to scale) of the 
TUNL-ITEP /3/3-decay apparatus. The enriched (98.4%) disk 
of 100 Mo is sandwiched between the two HPGe detectors. 

more or less matches the size of the two HPGe detectors. 
The disk thickness was chosen in order to achieve a suit- 
able balance between maximizing the source mass while 
limiting efficiency losses due to attenuation and geome- 
try effects. Preceding the study of this enriched disk, a 
metallic disk of natural molybdenum (1-kg mass, 10-cm 
diameter, 0.965-cm thickness, 9.6% 100 Mo) was studied 
to better understand potential backgrounds in the coin- 
cidence spectra. 

B. Efficiency Determination 

In order to properly interpret the 77 coincidence spec- 
tra obtained, a solid understanding of the efficiency of 
the apparatus to detect 77 coincidence events must be 
achieved. For the purpose of the present study, it was 
realized that a direct measurement of the efficiency was 
not only feasible, but it would also provide the most re- 
liable information. Using a source of known strength, 
the principle behind the efficiency measurement is to ob- 
serve a deexcitation cascade involving two 7 rays with 
energies and angular distributions comparable to those 
of interest from the /3/3-decay measurements. As such, 
the radioactive 102 Rh nuclide was judiciously selected for 
this measurement. 

The choice of the 102 Rh isotope is quite practical for 
many reasons. First, the half-life of 102 Rh is approxi- 
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FIG. 3: Decay scheme of Rh and the relevant gamma cas- 
cades of the excited states of 102 Ru. 



FIG. 4: (Color online) Calculated angular distribution, W , as 

a function of the angle 6 11 l2 between the two 7 rays emitted 

in the 77 cascade. The solid curve is the angular correlation 
for all + — > 2 + — > + decay sequences. The dashed (dotted) 
curve is the correlation between the two 7 rays emitted from 



the 21 



2+ -> 0+ sequence of 1U0 Ru ( 1U2 Ru). 



mately 207 days, which is convenient because the short- 
lived radioactivities decay away leaving a rather pure 
sample. Moreover, the source can be used over long pe- 
riods of time without having to apply significant correc- 
tions to a set of measurements performed over a period 
of a few weeks. Secondly, as shown in Figure|3 following 
the electron capture of 102 Rh to the 0+ state of 102 Ru 
with a branching ratio of 2.68%, the isotope emits two 
7 rays in coincidence with energies (468.6 keV and 475.1 
kcV) similar to those from the (3(3(0 + — > Of) decay of 
100 Mo. Since (3 + decay does not occur, there is no radi- 
ation from bremsstrahlung or annihilation, thus making 
the measurement very simple. At the same time, the 
probability to reach the Of state via electron capture 
to higher excited states is small (< 8.8%), hence most 
of these pairs of 7 rays (468.6 - 475.1 keV) are emitted 
without any additional quanta that could interfere with 
the measurement by a summation effect. 

In addition, the angular correlation of the two emit- 
ted 7 rays is identical to those of the 77 cascades that 
follow the (3(3 decay of 100 Mo to excited + final states be- 
cause it proceeds via the + — * 2+ — * + sequence. This 
condition arises because the corresponding deexcitation 
cascade proceeds solely through electric quadrupole (E2) 
radiation, with no mixing of other multipolarities. Fi- 
nally, one additional advantage of 102 Rh as the source is 
the presence of a higher excited 2 + state at 1103.2 keV 
of 102 Ru, corresponding to a subsequent 2 + — > 2 + — > + 
cascade that provides a measure of the detection effi- 
ciency for such a deexcitation sequence. Hence, the effi- 
ciency for detecting (3(3 decays to excited 2 + states can 
also be determined, albeit in a more complicated man- 
ner because the mixing of multipolarities (Ml and E2) 
depends on the specific 7 ray involved (see Figure QJ. 

A 102 Rh source was produced at TUNL by bombard- 
ing a target made from natural ruthenium (31.6% 102 Ru) 



with a 5-MeV proton beam, thereby generating 102 Rh via 
(p,n) activation of 102 Ru. Following the production of 
the 102 Rh source, a small portion of it (0.2cm x 0.2cm x 
0.1cm) was isolated for use in the efficiency measurement. 
Its strength was determined by comparing the intensity 
of the 475.1-keV gamma line to that of the 661.6-keV 
gamma line from a calibrated 137 Cs source. To repro- 
duce the photon attenuation in molybdenum that occurs 
in the /3/3-decay experiment, the source was placed be- 
tween thin natural molybdenum disks. Ten disks (10-cm 
diameter and 0.1-cm thickness) were used to ensure that 
the overall thickness was the same as the original exper- 
iment, thereby creating the same spacing between the 
two HPGe detectors. The use of these disks also enabled 
the variation of the position of the 102 Rh source with 
respect to the molybdenum disks, as a function of cylin- 
drical coordinates (radius, depth, and azimuthal angle) 
with respect to the center of the detectors. 

At each location of the 102 Rh source, the number of 77 
coincidence events involving the 469 - 475 keV cascade, 
arising from the decay of 102 Rh to the Of state of 102 Ru, 
were counted. The results were compared to the total 
number of decays that actually occurred, as given by the 
source strength and the branching ratio. Measurements 
of the 7 rays emitted from the 102 Rh source were adjusted 
to corresponding 100 Ru 7-decay sequences by applying 
corrections for both the relative efficiency of the HPGe 
detectors and the photon attenuation in Molybdenum. 
Figure [5] displays the results for the radial dependence 
of the efficiency £ 77 (r) to detect the 591 - 540 keV 77 
coincidence events from the 100 Mo (3(3 decay. [The nota- 
tion used henceforth, such as 591 - 540 keV coincidence, 
refers to 77 coincidence events in which a count in one 
detector is observed in the energy interval 591 ± 2.5 keV 
while at the same time a count is registered in the sec- 
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FIG. 5: (Color online) The efficiency e~ n (r) to detect 591 
- 540 keV coincidences, measured as a function of the radial 
distance r from the center of the detectors, and an appropriate 
fit to the data points. 



ond detector in the interval 540 ± 2.5 keV. The notation 
applies to all possible 77 combinations and energies, and 
the energy interval used (±2.5 keV) is reflective of the 
detector resolution.] 

An integration over the 100 Mo source volume yields the 
total detection efficiency and the results for the 77 
cascades of interest from excited states in 100 Ru are pre- 
sented in Table U for the enriched 100 Mo disk. A Monte 
Carlo simulation was also performed to check the validity 
of these values. It incorporated all known factors, includ- 
ing the full-energy peak efficiency of the HPGe detectors, 
the strongly anisotropic angular correlation of the 7 rays 
and their attenuation in the sample, and effects of the 
extended geometry. 



III. ANALYSIS 
A. Background Considerations 

Before proceeding with an analysis of the data from 
the TUNL-ITEP apparatus, it is important to address 
potential background contributions that may arise in the 
77 coincidence spectrum. Although the discussion that 
follows pertains to the (3(3 decay of 100 Mo to the 0^~ state 
of 100 Ru, the principle arguments can be easily extended 
to decays involving higher excited states (i.e., different 
7-ray energies). 

Due to the excellent energy resolution afforded by the 
HPGe detectors, the primary concern is background pro- 
cesses that can generate an identical gamma decay cas- 
cade with no additional quanta. Based on the present cri- 
teria to detect excited-state double beta decays by means 
of the simultaneous detection of the two 7 rays subse- 
quently emitted via deexcitation, any such background 
source would in principle be indistinguishable from the 



desired experimental signature. 

Realizing that the cascade involving two, and only two, 
7 rays with energies of 591 keV and 540 keV is unique to 
the deexcitation of the Of state of 100 Ru, the major back- 
ground consideration is of any mechanisms other than f3f3 
decay that may populate this excited level. Taking into 
account the fact that the 77 detection efficiency rapidly 
decreases as one moves away from the center of the 100 Mo 
disk, as seen in Figure 03 background sources that are of 
the greatest concern are those originating from within 
the source itself. 

As such, the most relevant background candidate 
would appear to be a significant flux through the detector 
setup of cosmic-induced protons, which can in principle 
result in the (p,n) activation of 100 Mo to bound states 
in 100 Tc. The subsequent (3 decay of this intermediate 
nucleus is delayed by the 15.8-s half-life, thereby ren- 
dering the anticoincidence counters ineffective in vetoing 
this process. The branching ratio of the 100 Tc (3 decay 
to the Of state of 100 Ru is 5.7%. Unfortunately, most 
of the remainder (93%) of 100 Tc (3 decays proceed to the 
ground state, which prevents using the 77 coincidence 
spectrum to search for additional gamma decay cascades 
that might otherwise provide a measure of the (p,n) ac- 
tivation of 100 Mo. 

An alternative means of estimating the proton flux 
in the apparatus and hence estimating the likelihood of 
background contributions from 100 Mo(p,n) is to search 
for the very similar (p,n) activation of 95 Mo to the isomer 
(1/2- , 61-day half-life) of 95 Tc. The electron capture of 
95m Tc populates two levels (1039.2 keV and 786.2 keV) 
of 95 Mo with large branching ratios (30.1% and 38.0%, 
respectively). The 1039.2-keV state decays to the ground 
state 89% of the time via a cascade of two 7 rays (835.1 
keV and 204.1 keV), while the 786.2-keV state proceeds 
to the ground state with a 77 cascade (where E 7 i = 582.1 
keV and E 72 = 204.1 keV) 78% of the time. The inter- 
mediate 204-keV level has a negligibly short half-life of 
0.75 ns. Searching for the 835 - 204 keV and 582 - 204 
keV coincidence events provides a convenient method for 
measuring the proton flux in the laboratory. 

To accomplish this study, a metallic disk of natu- 
ral molybdenum (1-kg mass, 10-cm diameter, 0.965-cm 
thickness, 15.9% 95 Mo) was investigated in the TUNL- 
ITEP apparatus for 180 days. Analysis of the data re- 
vealed exactly zero 835 - 204 keV coincidence events. For 
the 582 - 204 keV coincidence, there was no excess of 
counts observed beyond a fit to the background. Com- 
bining these two results allows one to gauge the potential 
for any contribution of (p,n) activation in the 100 Mo sam- 
ple, assuming a similar proton flux. The conclusion from 
the analysis of the 95 Mo study is that (p,n) activation 
contributes less than one count per year to the impor- 
tant 591 - 540 keV coincidence events of interest. 

Another potential background candidate involves 
sources arising from a significant neutron flux through 
the detectors and the 100 Mo disk. There exist two known 
mechanisms that could, in principle, provide unwanted 
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TABLE I: Calculations of the efficiencies £ l °^ (%) to detect 77 cascades resulting from the (3(3 decay of 100 Mo to various excited 
states of 100 Ru. 



Transition 



Level (keV) 



71 - 72 (keV) 



Branching Ratio fb 



-tot (0/ 



0+ 
0+ 
0+ 



or 
2+ 

o 2 + 



0+ - 0+ 



1130.3 
1362.2 
1741.0 

2051.5 

2387.4 



590.8 - 
822.4 - 
1201.4 

378.9 - 
1512.1 
689.4 - 
1847.8 
1025.2 - 



539.5 
539.5 

- 539.5 
1362.0 

- 539.5 
1362.0 

- 539.5 
1362.0 



1.00 
0.58 
0.59 
0.17 
0.86 
0.06 
0.51 
0.21 



0.219 ±0.022 
0.196 ±0.020 
0.185 ±0.019 
0.196 ±0.020 
0.165 ±0.017 
0.161 ±0.016 
0.142 ±0.014 
0.154 ±0.015 



77 coincidence events that would, in practice, be indis- 
tinguishable from the (3(3 decay to excited final states in 
the present experiment. 

First, if a ruthenium impurity resides in or near 
the 100 Mo disk, elastic scattering of neutrons on 100 Ru 
(12.6% abundance) could produce a 591 - 540 keV co- 
incidence that at present would register as a spurious 
100 Mo /?/3-decay event to the Of state of 100 Ru. A mea- 
sure of the significance of this possible background can 
be obtained by searching for similar neutron scattering 
processes, such as 100 Mo(n,n'7) and 76 Ge(n,n'j), where 
the abundance of 76 Ge is 7.44%. The former can result 
in an excited state of 100 Mo that decays to the ground 
state through a 77 cascade with E 7 i — 528.2 keV and 
E 7 2 — 535.6 keV, whereas the latter produces a 77 cas- 
cade with Eji = 545.5 keV and E~ (2 — 562.9 keV from 
the deexcitation of the 76 Ge excited state. 

During the experimental trial that studied the enriched 
100 Mo disk, approximately 12 events per year fit the 528 
- 536 keV coincidence profile. Also, for the neutron scat- 
tering process on the 76 Ge nucleus, the 546 - 563 keV 
coincidence was observed approximately four times per 
year. Making realistic assumptions about differences in 
cross section, detection efficiency, and a reasonable es- 
timate on an upper limit of the size of such an impu- 
rity within the 100 Mo disk, a conservative estimate of the 
threat of this neutron scattering background amounts to 
no more than approximately one half of one 591 - 540 
keV coincidence event per year. 

Secondly, there is a realistic, albeit unlikely, alternative 
scenario that could result in a coincidence involving 591- 
keV and 540-keV 7 rays that do not originate from the 
deexcitation of the 0^~ state of 100 Ru. Rather, following 
the neutron capture by a 100 Mo nucleus, 100 Mo(n,7), the 
(3 decay of 101 Mo can populate excited states of 101 Tc, 
from which two deexcitation gamma rays (Eji = 540.1 
keV and E l2 — 590.9 keV) are emitted as part of a four- 
photon cascade. 

In principle, summation and veto effects should rule 
out the likelihood of any significant contribution from 
this background candidate. Furthermore, the wide, 
neutron-broadened peaks that result from the neutron 
scattering interactions on germanium nuclei are not ac- 
companied by narrow peaks usually ascribed to thermal 



neutron capture. Thus, the thermal neutron flux is most 
likely very small. Even so, there is a small, non-zero prob- 
ability that the two bystander photons could escape un- 
detected, leaving behind a 591-keV photon in one HPGc 
detector and a 540-keV photon in the other. Conve- 
niently, an analysis of this possibility is readily accom- 
plished as a byproduct of the normal data compilation. 

The transition to the 2218-keV level of 101 Tc consti- 
tutes only 0.4% of 101 Mo (3 decays, and is followed 0.69% 
of the time by two possible four- photon cascades contain- 
ing both the 540.1-keV and the 590.9-keV 7 rays. A more 
likely scenario is the /3-decay to the 1319.6-keV level with 
a branching ratio of 6.8%, which produces a 77 cascade 
of two photons with energies of 713.0 keV and 590.9 keV. 
This 77 cascade is approximately 1000 times as likely to 
occur as the other four-photon cascades. An analysis of 
the 455 days of data from the 100 Mo run reveals four 
such 713 - 591 keV coincidence events. Thus, one can 
place a conservative upper limit of 0.003 counts per year 
originating from 540 - 591 keV background events from 
the four-photon deexcitation cascades involved in the (3 
decay of 101 Mo. Again, the actual contribution should, 
however, be considerably less because the two other pho- 
tons from the four-photon cascades would normally lead 
to the veto of the signal. 

Finally, for reference sake, it was found that truly ran- 
dom coincidences contribute less than 0.0003 counts yr _1 
for the 100 Mo experiment (with energies near 600 keV). 
Moreover, this background should be relatively uniform 
over the region of interest and would not be specific to 
the gamma decay cascades relevant to the present study. 



B. Decay to the 0^ Excited State 

The search for the (3(3 decay of 100 Mo to the 0^~ state 
of 100 Ru involves the simultaneous detection of the co- 
incident deexcitation (0^ — * 2f — > 0f s ) 7 rays, where 
E y i = 590.8 keV and E l2 = 539.5 keV (see Figure [TJ. 
This 77 cascade from the Of state of 100 Ru has a branch- 
ing ratio of 100%. The isotopically-enriched (98.4%) 
100 Mo disk was measured for a period of 455 days in 
the TUNL-ITEP apparatus. An inspection of the 7-ray 
spectra in coincidence with 539.5±2.5 keV and 590.8±2.5 
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FIG. 6: (Color online) The 7-ray spectra in coincidence with 
590.8 ± 2.5 keV (a) and with 539.5 ± 2.5 keV (b). The 22 
observed 591 - 540 keV events (shaded) from the /3(3(0 + — > 
Of) decay of 100 Mo were obtained in 455 days of measuring 
time. 

keV, respectively (as shown in Figure [BJ, reveals 22 such 
591 - 540 keV coincidence events. The width of the co- 
incidence window is five channels (5 keV) in order to 
account for the detector resolution. The presence of the 
peaks in their respective spectra is unmistakable when 
comparing their amplitudes to that of the surrounding 
background. A conservative measure of this background 
incorporates the full spectrum within ±100 keV of the 
591 - 540 keV coincidence peak, a measurement of which 
yields 2.50 ± 0.25 background counts per 5-keV bin. 

Deriving an empirical measure of the half-life of the 
excited-state (3(3 decay from the observed counting rate 
is straightforward. The 22 detected events (07^ = \/22 = 
4.7) yield iV 77 — 19.5 ± 4.7 events after subtracting 
the background, thereby corresponding to a signal-to- 
background ratio of approximately 8:1. From the sample 
size (1.05 kg), the counting time (~ 1.25 yr), and the de- 
tection efficiency (0.219%), the calculated decay half-life 
for the present measurement is 

T (0u+2u) = [ 6 . +i ; 9( s tat)±0.6(syst)] x 10 20 yr- (1) 

The asymmetrical statistical uncertainty results from 
the low number of counts, while the systematic er- 
ror entails the 10% uncertainty in the detection effi- 
ciency. It should be noted that this value for the half- 
life differs slightly from the previously published value of 



[5.9+^(stat) ± 0.6(syst)] x 10 20 years 23] for this exper- 
iment. The difference arises from a slight modification to 
the calculated efficiency, the inclusion of some additional 
data (455 days of counting versus the original 440 days), 
and a corrected calculation of the statistical error bars. 

From the decay half-life presented here, it is possible 
to extract a value for the nuclear matrix element corre- 
sponding to this 2v(3(3 transition. Using the phase-space 
factor G = 1.64x 10 -19 yr -1 (for the axial-factor coupling 
constant gA — 1.254), one obtains the nuclear matrix el- 
ement M2v{$i) = 0.101 ±0.013, as scaled by the electron 
mass. 

The detection scheme is incapable of deciphering the 
two-neutrino mode from the neutrinoless mode of (3(3 de- 
cay, thus the half-life of Equation ^ is listed as an in- 
clusive {Qv + 2v) quantity. Nevertheless, using values 
for 100 Mo (3(3 transitions to the ground state of 100 Ru 

( T i/2 = 6 - 75 x 1018 y r and T i/ 2 > 5 - 5 x 1022 

yr, 90% C.L. |2{|), one concludes that the neutrinoless 
mode constitutes less than 0.012% of (3(3 transitions to 
the ground state. In fact, the most recent 100 Mo data 
from the NEM03 experiment [2(j restrict the neutri- 
noless (3(3 decay of 100 Mo to the ground state further: 
T i/2 = I™ ± 0.02(stat) ± 0.54(syst)] x 10 18 yr and 
Tjy 2 > 3.1 x 10 23 yr (90% C.L.). Thus, assuming a similar 
proportion for the excited-state (3(3 decays, it is reason- 
able to assume a negligible contribution from 0v events 
in the present data. This conclusion is further reinforced 
by theoretical predictions for the neutrinoless (3(3 decay 
of 100 Mo to the 0i excited state, with the half-life for this 
0v(3(3(p + -> 0^) transition estimated to be: 2.59 x 10 26 
yr [23 and (0.76 - 1.46) x 10 25 yr for (m„) = 1 eV. 

In Figure \7\ the 7-ray spectra in coincidence with 
539.5 ± 2.5 keV and 590.8 ± 2.5 keV show a broader 
range of energies and reveal some interesting features 
of the background. Both spectra exhibit the presence 
of some known gamma peaks, such as the annihilation 
peak (511 keV) and peaks from 208 T1 (583 keV) and 214 Bi 
(609 keV). As expected, the background at lower ener- 
gies (< 500 keV) is generally more prominent. This phe- 
nomenon is a consequence of two factors: (1) the larger 
photopeak detection efficiency at lower energies and (2) 
the contributions at low energies from Compton scatter- 
ing events. In the latter case, the Compton continua of 
a large number of the more prominent primordial back- 
ground 7 rays, in addition to those of interest in the 
present study, lie below 500 keV. 

Another interesting feature that was observed in the 
77 coincidence spectra at higher energies (> 800 keV) is 
the Compton scattering of at least two known primor- 
dial background peaks. The most noticeable instance 
involved the 1461-keV 7 ray ( 40 K), which is always emit- 
ted in single with no accompanying gamma quanta. The 
Compton scattering of this photon can lead to the depo- 
sition of 591 (540) keV in one of the HPGe detectors and 
the remaining 870 (921) keV in the other HPGe detector, 
and there are indeed rather prominent peaks at the re- 
spective energies of the corresponding spectra of Figure 
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FIG. 7: (Color online) The 100 Mo 77 coincidence data, compressed to include 5 keV per bin. Shown are the 7-ray spectra in 
coincidence with 590.8 ± 2.5 keV (a) and with 539.5 ± 2.5 keV (b), with some of the prominent background lines identified (CS 
= Compton Scattering). Note the 22 observed 591 - 540 keV coincidence events from the (3f3(0 + — > Of) transition. 



The same phenomenon occurs for the 1764.5-keV 7 
ray ( 214 Bi), producing the observed 1225 - 540 keV and 
1174 - 591 keV coincidence events. To a lesser extent, 
there also appear to be some Compton scattering coinci- 
dence events corresponding to 1001.0-keV ( 234m Pa) and 
1120.3-keV ( 214 Bi) 7 rays. 



C. Decay to Higher Excited States 

In addition to the transition to the excited 0^~ state, 
the (3(3 decay of 100 Mo can also populate higher excited 
states of 100 Ru. In principle, the measurement of these 
additional transitions is readily accomplished with the 
TUNL-ITEP apparatus by searching for alternative 77 
coincidence events associated with the deexcitation of 
these higher excited states. In the present study, the 
relevant 100 Ru levels of interest and the corresponding 
gamma decay cascades are shown in Figure |H1 

Since these other processes were not observed in the 
present experiment, it is useful to briefly mention the 
procedure for setting a lower limit on a decay half-life. 
In general, the half-life limit is given by 



T1/2 



> 



(\n2)N tf b e tot 



77 



N d 



(2) 



where ft denotes the branching ratio of a given 77 cas- 
cade, and all relevant values for the detection efficiency 
are taken from Table||| As recommended by the Par- 
ticle Data Group |29j, the factor N4, which represents an 




(keV) 



FIG. 8: Decay scheme for the (3/3 decay of 100 Mo to higher 
excited states of 100 Ru, excluding the 0^ state. Shown are 
the subsequent gamma cascades and the corresponding 7-ray 
energies and branching ratios. 



upper limit on the number of detected events above back- 
ground, is determined by a statistical estimator derived 
for a process that obeys Poisson statistics. 

The deexcitation of the higher excited + states of 
100 Ru at 1741.0, 2051.7, and 2387.2 keV can each oc- 
cur via two different 77 cascades, one involving the tran- 
sition through the 2^ state (539.5 keV) and the other 
through the 2 J state (1362.0 keV). Hence, the search 
for (3(3 decay to these excited states involves monitoring 
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the 7-ray spectrum in coincidence with these two 7-ray 
energies. It should be noted that, in addition to these 
two-photon cascades, there is also a three-photon deexci- 
tation sequence. However, the existing apparatus is not 
sensitive to cascades involving more than two 7 rays, thus 
such decay sequences will not be considered here. At the 
same time, a search for the (3(3 decay to the 2^~ excited 
state is accomplished by searching for coincidence events 
involving 822.4-keV and 539.5-keV 7 rays, without the 
emission of any additional gamma quanta. 

An analysis of the 455 days worth of data yielded the 
following results for the (3(3 decay of 100 Mo to higher ex- 
cited states. For the transition to the 0^" state (1741.0 
keV) of 100 Ru, zero 1201 - 540 keV coincidence events 
were detected and one 379 - 1362 keV event was observed, 
which corresponds to a counting rate that is consistent 
with the background. In the case of the transition to 
the 0% state (2051.7 keV), a search revealed one 1512 - 
540 keV event (again consistent with the corresponding 
background rate) and zero 689 - 1362 keV events. For 
the transition to the Oj state (2387.2 keV), zero 1848 
- 540 keV events and zero 1025 - 1362 keV events were 
observed. Lastly, a search for 822 - 540 keV coincidence 
events, corresponding to the (3(3 decay to the 2^~ state, 
yielded zero counts. Taking into consideration the associ- 
ated 77 background rate in the regions of interest, limits 
on the /3/3-decay transitions to higher excited states were 
calculated from Equation [3 and are listed in Table HTl 



IV. CONCLUSIONS 

The implementation of a 77 coincidence counting tech- 
nique has resulted in a half-life of the (3(3(0 + — > 0+) tran- 
sition of 100 Mo of [6.0ti;?(stat)±0.6(syst)] x 10 20 yr, 
and the corresponding nuclear matrix element is calcu- 
lated to be M 2v (0t) = 0.101±0.013, as scaled by the elec- 
tron mass. This half-life is in very good agreement with 
previous results for this transition [20ll2l| . and represents 
an independent confirmation of this decay. Moreover, the 
success of the present experiment and the methods used 
is exemplified by a superior signal-to-background ratio. 
Although the decay to additional excited final states was 
not observed in the present measurement, stricter lim- 
its on the decay half-lives for these transitions were ex- 
tracted. 

This result is particularly important in the context of 
the single state dominance hypothesis (SSDH) [33, El 
133 Using the result for the EC transition 100 Tc 100 Ru 
[33, a prediction for the (3(3 decay of 100 Mo to the 0^ 
state of 100 Ru has been made: if" = 4.45 x 10 20 yr 
|34|. Unfortunately, this prediction has an accuracy of 
only 50% due to the uncertainty of the measured EC 
transition. However, if in the future there is significant 
improvement in the accuracy of the EC transition, it may 
be interesting to compare the experimental half-life as 
presented here to the theoretical prediction to check on 
the validity of the SSDH. 



On a more general note, it is interesting to point out 
that the present 77 counting technique represents a pow- 
erful tool that could potentially be extended to future, 
large-scale neutrinoless double beta decay experiments, 
such as Majorana [H, GENIUS HH, SuperNEMO 
and CUORE [jj^, for example. Although Qv(3(3 transi- 
tions to excited final states have smaller values, the 
emission of fixed-energy, deexcitation 7 rays opens the 
possibility of detecting the 0v(3(3 decay by way of a mul- 
tiple, unambiguous coincidence involving the successful 
detection of both the emitted (3 particles as well as the 7 
rays. Such a multiple coincidence technique would result 
in a dramatically-reduced background, provided the de- 
cay products are detected with high efficiency and with 
good energy resolution. 

Indeed, the information that can be gained from 
searching for neutrinoless double beta decay to excited 
final states is just as valuable as that from the transi- 
tion to the ground state. Furthermore, depending on the 
particular experimental design and apparatus, it may in 
principle be possible to incorporate the option to search 
for Qv(3(3 transitions to excited final states, in conjunc- 
tion with the search for the transition to the ground state, 
with minimal cost or reconfiguration. Adding this pos- 
sibility would be considerably advantageous as it would 
serve to greatly augment the breadth of a given experi- 
ment, thereby increasing the potential for discovery. 
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TABLE II: A summary of the experimental half-life values, or limits, for the f3(3 decay of 100 Mo to various excited states of 
100 Ru corresponding to a 455-day counting period. Limits are indicated at the 68(90)% C.L. References for previous results 
are given in square brackets. 

Transition Level (keV) Q 0I3 (keV) t[°" 2 +2v) (xlO 20 yr) Previous results (xlO 20 yr) 

0+ 0+ 1130.3 1904.0 6.0l^(stat) ±0.6(syst) 6.1+_^ [20] 

9.3ti : ? [21] 

0+ -> 2+ 1362.2 1672.1 > 54(27) > (13) [20] 

0+ -> Oj 1741.0 1293.3 > 55(28) > (13) [20] 

0+ -> 0t 2051.7 982.6 > 42(24) 

0+ -> 0t 2387.2 647.1 > 49(25) 
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